Lateral migration of a 2D vesicle in unbounded Poiseuille flow 
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The migration of a suspended vesicle in an unbounded Poiseuille flow is investigated numerically 
in the low Reynolds number limit. We consider the situation without viscosity contrast between 
the interior of the vesicle and the exterior. Using the boundary integral method we solve the corre- 
sponding hydrodynamic flow equations and track explicitly the vesicle dynamics in two dimensions. 
We find that the interplay between the nonlinear character of the Poiseuille flow and the vesicle 
deformation causes a cross-streamline migration of vesicles towards the center of the Poiseuille flow. 
This is in a marked contrast with a result [L.G. Leal, Ann. Rev. Fluid Mech. 12, 435 (1980)] 
according to which the droplet moves away from the center (provided there is no viscosity contrast 
between the internal and the external fluids). The migration velocity is found to increase with the 
local capillary number (defined by the time scale of the vesicle relaxation towards its equilibrium 
shape times the local shear rate) , but reaches a plateau above a certain value of the capillary num- 
ber. This plateau value increases with the curvature of the parabolic flow profile. We present scaling 
laws for the migration velocity. 
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I. INTRODUCTION 

Vesicles are closed phospholipid membranes suspended in an aqueous solution. They constitute a first step in a 
model aiming to capture elementary ingredients in the study of the dynamics of more complex entities such as red 
blood cells. Of particular interest is the migration of blood cells in the circulatory system. The study of migration 
of soft particles under flow presents fundamental (e.g. understanding the subtle interplay between deformation and 
the flow) as well as technological interests (e.g. understanding this problem may help monitoring vesicles and cell 
guidance in various circumstances, such as in microfluidic devices, in the process of cell sorting-out, and so on...). 
In the present work we focus our attention on describing the dynamics of a single suspended vesicle in a non-linear 
shear gradient of a plane Poiseuille flow. We consider the small Reynolds number limit, so that inertia can be 
neglected. Vesicles in flow field have been the subject of extensive studies, both in an unbounded linear shear flow 
P, 13, H, S, H fl S [3 ^ ^^11 ^ the presence of a wall [^, [l^, [HI ■ In an unbounded linear shear flow (of low 
Reynolds number) a vesicle does not exhibit a lateral migration with respect to the flow direction. The presence of a 
wall breaks the translational symmetry perpendicular to the flow direction and a vesicle is found to migrate away from 
the wall [1, S [O, ■ This so called lift force is caused by the flow induced upstream-downstream asymmetry of 
the vesicle Q. More recently, it has been shown that even a spherical vesicle may execute a lift force as well, provided 
that the wall is flexible [l^. In that case, the wall deformability breaks the upstream-downstream symmetry. 

A nonlinear shear flow has a non translationally invariant shear rate. It is therefore of great interest to understand 
its possible contribution to a cross-streamline migration process. Here we focus on the pure bulk effect due to the 
non-linear shear gradient of a plane Poiseuille flow alone. Accordingly, we consider a parabolic flow profile in the 
absence of bounding walls (in order to avoid any wall induced lift force) , and we consider neutrally buoyant vesicles 
so that gravity effect is suppressed. 

We find that the curvature of the imposed velocity profile, together with the vesicle deformability, causes a systematic 
migration of a tank-treading vesicle perpendicularly to the parallel streamlines towards the flow center-line. Our results 
show that the migration velocity increases with the curvature of the flow profile and we provide also scaling results for 
the migration velocity as a function of the local capillary number. This behavior is different from a prediction made 
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for droplets according to which droplets should migrate away from the center of the Poiseuille flow towards the 
periphery. Actually, in Ref. UJl] it is predicted that the direction of the lateral migration of a droplet depends on the 
viscosity contrast. For values between 0.5 and 10 -and particularly in the absence of a viscosity contrast as treated 
here-, migration occurs towards the periphery, while for values smaller than 0.5 or greater than 10, it occurs towards 
the center line. We did not observe any of these scenarios neither from numerical studies. 

The scheme of this paper is as follows. In Sec. 2 we present the model equations and describe briefly the method 
used to solve the problem. In section 3 we define the dimensionless parameters and provide typical experimental 
values. In Section 4 numerical results and their discussion are presented. Section 5 is devoted to some concluding 
remarks. 



II. MODEL AND METHOD 



A. Hydrodynamical equations and boundary integral method 

The flow of an incompressible Newtonian fluid with viscosity 77 and density p is characterized by the dimensionless 
Reynolds number, 

Re=^, (1) 
V 

where Vq is a characteristic velocity and Rq a characteristic length of the studied system. In our case we take the size 
of a vesicle, which is of the order of 10 — lOO/im T^, as the characteristic length. For such a length and for vesicles 
suspended in an aqueous solution subject to shear, with moderate applied shear rates (7 = Vq/Ro) that are usually 
of the order of 10s~^, the Reynolds number is rather small. Re ~ 10^^ — 10^"^ ^ 1. Therefore, the flows of the 
fluids inside and outside the vesicle, which are taken to be of the same nature, are well approximated by the Stokes 
equations, 

-Vp + ?7V2v = f, 
V • V ^ 0, 

where v is the fluid velocity, p is the pressure, and f is the force imposed by the deformable vesicle membrane on the 
two fluids (it is a local force having a non zero value only at the membrane). This force is given by the functional 
derivative of the vesicle membrane energy with respect to the membrane displacement, as discussed in the next 
subsection. 

Using the boundary integral method [lBili3| we solve the equations ^ in two-dimensional space. The velocity field 
at any point in the fluid (or at the membrane; the membrane velocity is equal to that of the adjacent fluids provided 
that the membrane is not permeable, and that there is no slip at the membrane) can be written as a superposition 
(due to linearity of the Stokes equations) of two terms, namely the contribution from the vesicle boundary, plus a 
contribution due to the undisturbed applied Poiseuille flow v(r)pois (to be written below), 

v,{r)= f dr'Gy(r-r')/j(r')ds(r')+t'^(r)pois. (3) 
Jan 

Here dfl refers to the vesicle boundary. Gy denotes the Oseen tensor, also called Green's function of the Stokes 
equations. Since we focus on dynamics of a vesicle suspended in an unbounded domain, we use the two-dimensional 
free space Green's function, that has the following expression [l3, HO] : 



G.=^(-%ln. + ^), (4) 

where r = |r — r'| and is the i*'' component of the vector r — r'. 

Equation ^ is valid in the fluid as well as at the membrane. In order to obtain the membrane velocity we replace 
r by the membrane vector position. Numerically, the vesicle membrane contour (in 2D) is discretized, as explained 
in Ref. [l^l . After evaluating the membrane force which enters the right hand side of Eq. ([3]) (see next section for 
the force evaluation), the velocity is then evaluated at each discretization point using Eq. ([3]). The displacement in 
the course of time of the vesicle membrane is obtained by updating the discretization points after each time iteration, 
using an Euler scheme, r{t + dt) = v(r, t)dt + r(t). In the following section we shall discuss in more detail the forces 
and the constraints. 
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B. Membrane force and comparison with droplets 

The vesicle membrane is a bilayer made of phospholipid molecules having a hydrophilic head and two hydrophobic 
tails. At room temperature (and at physiological temperature as well) the membrane is fluid. The membrane can be 
viewed as a two-dimensional incompressible fluid. This incompressibility property implies the inextensibility of the 
membrane, and therefore, the conservation of local area. Moreover, since the vesicle encloses an incompressible fluid 
and the membrane permeability is very small, the vesicle volume must be a conserved quantity. Due to membrane 
impermeability, the membrane velocity is equal to the fluid velocity of the adjacent layer. Consequently, and because 
we use explicitly the incompressibility condition V.v — for fluids, the enclosed volume is automatically conserved. 
The area of the membrane is not conserved automatically (think of a droplet that can spread out on a substrate; 
its volume is conserved while its surface increases). Thus, in order to fulfill membrane local area (or perimeter in 
2D) conservation, one must introduce a surface (local) Lagrange multiplier. This Lagrange multiplier C(s, t) depends 
on the curvilinear coordinate s along the vesicle contour and on time, since the incompressibility should be fulfilled 
locally and at each time. C is the surface analog of the pressure field p(r, t) which enforces local volume conservation 
of a 3D fluid. 

Due to the fact that phospholipidic molecules are bound to the membrane (there is no exchange between the bilayer 
and the surrounding solution), the local area (or area per molecule) remains constant in the course of time. This is 
a major difference with droplets, where the surface molecules can easily migrate into the underlying bulk and vice 
versa, causing thus a change of area. One may say that the surface molecules of a droplet are in contact with a large 
reservoir of molecules in the bulk (the reservoir fixes the chemical potential, while the number of molecules at the 
surface is a fluctuating quantity). For a droplet, the question of energy per unit surface makes sense: it corresponds 
to the increase in energy due to the transfer of a molecule from the bulk towards the surface. This is surface energy, 
or surface tension (note that for a liquid the surface energy and surface tension refer to the same quantity, while this 
is not the case for a solid). The surface force for a droplet is given by the classical Laplace law 

id = -aHn (5) 

where a is the surface tension, n is the outward normal unit vector, and H is the surface mean curvature. Note, that 
in the 2D case, as treated here, there is only one curvature H and H is by convention counted to be positive for a 
circle. 

As a vesicle membrane does not naturally change its area, the notion of cost of energy per unit area can not be 
evoked. Moreover, one may think of changing the area by applying an external force and the surface energy is in 
this case the work associated with the applied force. This notion is, however, quite different from surface energy of 
a droplet. Here, for vesicles, we must rather refer to a surface stress, since by applying a force the intermolecular 
distance (due to stretching, for example) is modified. In contrast, a droplet may change its area without affecting the 
intermolecular distance. If we do not apply a large force in order to stretch the microstructure of the lipid layer, then 
membrane incompressibility is fulfilled. Hydrodynamical forces (as those encountered in this problem) are too small 
in comparison with cohesive forces, so that membrane incompressibility is safely satisfied. 

From the mechanical point of view, the membrane can be viewed as a thin plate, where the soft (or easy) mode is 
the bending one. The corresponding energy is given by the Helfrich curvature energy This reads in 2D 

Ec = ^ f H^ds, (6) 
^ Jon 

where k is the membrane rigidity and H the local membrane curvature, ds is the elementary arclength along the 
vesicle contour. Note that for the sake of simplicity, we do not account for a spontaneous curvature (a constant 
spontaneous curvature, Hq, may be included by substituting H hy H — Hq). In order to take into account the area 
(perimeter in 2D) constraint, we must add to the above energy the following contribution J^^^ ^(s, i)ds, so that the 
total energy reads 

E^^ f H^ds+ f as,t)ds (7) 
^ Jon Jdn 

where C,{s,t) is a local Lagrange multiplier. Note, that global conservation of the perimeter would be unphysical, 
because it would allow at some range of the membrane an arbitrarily large stretching and at the same time at another 
point a corresponding compression in a way that the global length is preserved. As discussed above, stretching or 
compression is possible only under the action of strong forces, of the order of cohesive forces. 

The force acting on the membrane is obtained from the functional derivative of the vesicle energy E with respect 
to a membrane displacement. The resulting force has been already used previously (e.g. ^8, _20]), but a derivation has 
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not been presented. A detailed derivation is given in the appendix. The resulting force is 



where t is the tangent vector (and recall that n) is the normal vector). This force is composed of a normal as well 
as a tangential contribution. If C is constant, then only the normal part survives because of the following reason. 
If C is constant, the tension-like force (which is a vector) associated with C, is tangential to the curve, and has the 
same magnitude at both extremities of an arc element ds (which can be taken to be a portion of a circle, provided 
that ds is small enough). It follows, that the sum of the two forces is directed in the normal direction. If, on the 
contrary, C changes along the contour, then the two values at the extremities of ds arc different, and the force has, 
besides a normal part, a tangential one, which is given by {dC/ds)t. On the other hand, the bending energy depends 
on the curvature (which is a geometrical quantity). It follows that the only force that is able to change the shape of 
a geometrical surface (i.e. a mathematical boundary having no internal physical structure) must be normal to the 
surface. Finally, note that the term —(Hn has the same structure as the force due to surface tension of a droplet 
Eq. ([5]). There is, however, a significant physical difference: for a droplet a is an intrinsic quantity which represents 
the cost in energy for moving a molecule from the bulk (surrounded by other molecules) to the surface (and thus it 
looses some neighbors). In the present problem C is a Lagrange multiplier which must be determined a posteriori by 
requiring a constant local area. C is not an intrinsic quantity, but rather it depends on other parameters (like n, the 
vesicle radius, etc.). 



C. Fulfilling local membrane area 

In principle, from Eq.([3]) we can determine the membrane velocity, if the force and the initial shape are given. The 
force (|S1) contains geometrical quantities (like the normal and H) which are determined from the initial shape, plus 
a function ({s,t), which is unknown a priori. Numerically, the following method has been tested. An initial shape 
(typical an ellipse) and an initial (typical constant along the contour) have been chosen. Then the geometrical 
quantities appearing in the force can be calculated (the method of discretization of the integral equation ^ has been 
discussed in [20|). This allows one to evaluate the right hand side of ^ at initial time. The membrane velocity at 
this time is thus fixed. We then displace each membrane element according to the computed velocity, and by this way 
we obtain a new shape. However, the new shape does not fulfill, in general, the local membrane incompressibility. 
A local stretching (or compression) of the membrane takes place as long as the projected divergence of the velocity 
field of the fluid adjacent to the membrane is non zero. We thus must adjust the appropriate function (^(s) in order 
to fulfill this condition. The condition that the projected divergence must vanish reads 

(I - nn) : Vv = (9) 

where I is the identity tensor, and nn stands for the tensor product (I — nn is the projector on the contour). The 
above relation can be viewed as an implicit equation for C{s), similar to V.v = which fixes the pressure field in 3D 
fluids. This way of reasoning is quite practical in the analytical study of vesicles [lij . From the numerical point of 
view, this way has suffered from several numerical instabilities. We have thus introduced another approach [20| as 
outlined below. 

In a 2D simulation, when discrctizing the vesicle membrane contour, the vesicle perimeter conservation constraint 
could be achieved without dealing with the local Lagrange multiplier entering the membrane force given by Eq. ([H]). 
This constraint could be fulfilled in another and more convenient way. For that purpose we have used a straightforward 
method based on the fact that two material representative points on the membrane are attached to each other by 
strong cohesive forces which we describe by quasi-rigid springs, so that we can achieve in numerical studies less than 
1% variation of the area. By this way an additional parameter kg is introduced, which is the spring constant [20j . 
C^(*) = ks{ds{i) — (is°(i)). By choosing kg large enough (in order to keep the membrane quasi-incompressible) the 
discretization step ds{i) is kept as close as possible to its initial value ds^{i). Typically in units where rj = k = 1 and 
where the typical radius of vesicles is of order unity, a value of ks = 10^ has proven to be sufficient. 



D. Applied flow 



The applied Poiseuille flow v(r)pois has the following form 



^x('^)pois — ^max 
Wy(r)pois = 0, 



1 



fy_Y 



(10) 



5 



where t^max is the maximum velocity at the centerhne located at y = and 2'w is the width of the Poiseuille profile. 
For our simulations we choose always the aspect ratio Rq/w 1 in order to keep v-^{w) = v-^i—w) — practically 
unperturbed by the presence of the vesicle. 



III. DIMENSIONLESS NUMBERS 



It is convenient to use in the simulation a dimensionless parameter that we call the local capillary number, which 
we define as 

C,(r)=r7(r). (11) 

r is the characteristic time for a vesicle to relax to its equilibrium shape (in the absence of imposed flow) , which is 
given by 

r = ^, (12) 

K 

7(r) is the local shear rate of the applied Poiseuille flow, that can be evaluated from the corresponding velocity profile 

Here c is the curvature of the Poiseuille flow profile, which is given by c = d^Vx/d'^y. In the numerical scheme, there 
is another capillary number associated with the tension fcg (or spring constant), and is defined by Cas = TjRo/ks. In 
most simulations we have kept the ratio Ca/Cas small (of the order of 10"'^). This means that the time scale for 
stretching/compression of the membrane is fast in comparison to bending. In other words, local area conservation is 
adiabatically slaved to the overall shape evolution. 
As a characteristic velocity we choose. 



^o = ^ = r^, (14) 



with i?o = L/2'K, where L is the vesicle perimeter. Hereafter we shall use r as a unit of time, i?o of length and 
Vq as unit of velocity. For typical experimental values of rj (e.g. water), and by using standard values for vesicles 
K, ~ 20kBT {ksT is the elementary thermal excitation energy) and Rq ~ lO/im one finds r 10s and Vq ~ 1/xm/s. In 
the following (and especially in the figures of the simulation), when a velocity is written in terms of a number (without 
units) this means it is expressed in unit of Vq- Since Vq is typically of order lfj.m/s, the velocity is given practi cally 
in iim/s. The reported values for the velocities in experiments on vesicles are in the range of 1 — lOOfim/s (gI. |24|. 
In the circulatory system, data are for the shear stress at the vessel wall are well documented [2l|. For example, in 
arteriesf2lj the shear stress is of about 1 — 2 Pa. Dividing this by the Plasma viscosity (close to that of water), one 
finds the shear rate at the wall, jwaii ~ 10^s~^. The velocity at the center of the arteries is of about 'Ywaiiw. For small 
arteries w ~ 100 /im, so that Vmax ^ 10^ fj,m/s (for venules, one has about Vmax ~ 10'*/im/s). The chosen values in 
the simulations (see figures in the next section) are rather in the experimental range for vesicles, but are not far away 
from data on blood flow in arteries. 



IV. SIMULATION RESULTS AND DISCUSSION 



The left part of Fig. [T] illustrates a free vesicle in an unbounded plane Poiseuille flow, the dynamics of which is 
investigated. The right part of Fig. [T] shows the time evolution of the lateral position of a vesicle which has been 
released initially from five different vertical positions: yo = 0, 1,2,3,4. In most cases we have studied quasi-circular 
vesicles placed in a Poiseuille fiow characterized by Wmax = 800 and w = 10. In all situations treated so far, vesicles 
migrate to the center of the Poiseuille flow where no further lateral migration is observed. The position gives the 
distance from the center of the Poiseuille flow measured in units of the vesicle effective radius (as defined in the last 
section). All the curves are linear in a large range and deviations from this linear law occur only close to the center 
of the Poiseuille flow, in a range smaller than the vesicle size. 

During the migration the vesicle shape undergoes deformations due to the hydrodynamic stresses imposed by the 
Poiseuille flow on the membrane. The vesicle is deformed and tilted until reaching a quasistationary orientation which 
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FIG. 1: (Color) The left plot shows a schematic representation of the migration of a vesicle in an parabolic flow proflle 
corresponding to an unbounded plane Poiseuille flow. The right plot shows the time evolution of the lateral position of a vesicle 
released initially at flve different initial positions yo = 0, 1, 2, 3, 4 in units of the characteristic time r (right). 




FIG. 2: (Color) The shape of the vesicle changes from an initially elliptical shape in part a) to the final parachute shape in d) 
when it migrates towards Poiseuille flow center-line, reduced area v — 0.90. 

is oblique with respect to the parallel streamlines. Figure [2] shows different vesicle shapes deformation occurring in 
our simulations during the migration, from an initially elliptical shape at the initial position yo = 1 shown in Fig.[2l[a), 
to a final parachute shape at the center of the Poiseuille flow as shown in Fig. [21[d). More or less similar parachute 
shapes are known for capsules and red blood cells [23, l23| as they have been observed also experimentally for vesicles 
in Ref. ^23|, but all these examples concern capillary flows. Thus, the present result shows that the parachute is not 
necessarily enforced by a wall, but rather by the curvature of the Poiseuille flow itself. 

Before the vesicle reaches the center of the Poiseuille flow, it acquires an asymmetric shape as depicted in Fig. [Hb) 
and in Fig. [2lc). This asymmetry, which is caused by the non- uniform shear rate across the vesicle, is crucial for 
cross-stream line migration of vesicles in a plane Poiseuille flow. 

As stated above, a drop (having no viscosity contrast with the ambient fluid) is predicted to drift towards the 
periphery [l3 |. Thus, vesicles and droplets behave quite differently. In section IIB we have presented the main 
differences between vesicles and droplets, both from the physical and the mathematical point of view. Which of these 
differences, albeit very important, may explain the differences in the migration direction, is not clear at present. For 
vesicles, we have explored a large domain of parameter space and in all cases the vesicle migrate towards the center. 

The migration velocity depends on various parameters. Of particular importance are the curvature of the velocity 
proflle of the Poiseuille flow and the local capillary number, as discussed in the following. To the best of our knowledge, 
there is in the absence of a wall no lateral migration in a linear shear flow. In the presence of a flow !!! with a 
nonlinear shear gradient, migration becomes possible, provided that the shear rate changes on the scale of the vesicle 
size. Therefore, curvature of the Poiseuille flow proflle plays an essential role, but more precisely, the magnitude of 
the local capillary number, which determines essentially the vesicle deformation (which loses the up-down symmetry 
due to the shear gradient), is the most relevant quantity. 

The dependence of the migration velocity on the local capillary is shown in Fig. [3] for different values of fmax, w 
and c, after the decay of an initial transient. In Fig. [3l^a) and in Fig. [3jb) we kept the value of Umax fixed and we 
investigated the vesicle migration by varying the value of w. For smaller values of w, which corresponds to larger 
values of the curvature c, the vesicle migrates faster towards the center of the Poiseuille flow. Fig. ^h) shows the 
data collapse by plotting the migration velocity normalized to the curvature versus the local capillary number. In 
Fig. [3jc) and in Fig. [3fd) we kept w fixed and we examined the effect of varying the value of i^max for each value 
of w. The vesicle migrates faster with increasing values of v^ax for every fixed w, because the curvature c increases 



7 



a) 




= 500 




w - 


5, c = 


40.00 




w = 


6, c = 


27.78- 




w - 


7, c = 


20.40 . 













01234567S9 10 
Time 



0.020 - 


b) 




0.015 - 






0.010 - 
0.005 - 




v^, = 500 

w = 5. c = -40.00 

w = G.c = -27.78 

w=7, c= -20.40 


0.000 - 







20 40 60 80 100 120 
Local Capillary Number 





20 40 60 80 100 120 
Local Capillary Number 



e) 


1 1 1 1 1 1 1 1 1 1 1 1 1 

c =-40 




v^^ = 500, w = 5 ■ 




v^^ = 720, w = 6- 




= 9S0, w=7 ' 









3 4 5 6 
Time 




20 40 60 80 100 120 
Local Capillary Number 



FIG. 3: (Color) Time evolution of the vesicle position in an unbounded Poiseuille flow and its corresponding normalized 
migration velocity versus the local capillary number for different values of Wmax, w and c (see text). The data correspond to 
situation where initial transients have decayed. 



with Wmax- Data collapse is again obtained in Fig.[3l^d) by plotting the normalized migration velocity versus the local 
capillary number. The data collapse is more pronounced for smaller values of the curvature. In Fig.[3je) and Fig. [3ljf) 
we have varied Wmax and w in such a way to keep the curvature fixed. We find that the vesicle migrates in this case 
to the Poiseuille flow center-line for the three parameters combinations exactly (i.e. quantitatively the same results) 
in the same manner, which emphasizes again the important role of the nonlinear shear fleld. From the above study, 
we can conclude that the migration velocity in an unbounded Poiseuille flow normalized to the curvature c should be 
described by the following universal scaling law: 

''~"°"^^^ -/[g.(y)]. (15) 
c 

The extraction of this law is based on results of Figs. [Sb, EH, and[3f. The function / is universal and depends only 
on Ca- The analytical form of the universal function is not at present. A first step towards this issue is to develop an 
analytical theory in the small deformation limit, as in Ref.[l9'|. The small deformation limit provides us with nonlinear 
differential equations for the shapes and the migration of the vesicle instead of the less tractable integro-differential 
equation ([3]). With this semi-analytical approach progress seems more likely for both, for an understanding of the 
migration direction, and the determination of the scaling function /. 

If the initial vesicle shape is not quasi-circular but elliptical, we find a similar behavior as depicted in Fig. [31 
The deformability of the vesicle, which depends on the bending rigidity k, is a further ingredient for migration. By 
increasing the bending rigidity k the local capillary number Ca oc decreases and so does the migration velocity. 
This leads also to the conclusion that a rigid particle, corresponding to very large values of k, will not exhibit a lateral 
migration in parabolic shear flow in the Stokes limit. It has been shown earlier that rigid spheres migrate only due to 
the contribution of {v ■ V)f in the Navier-Stokes equation 1251 . which is beyond the Stokes limit. Similar trends as for 
a vesicle are obtained for deformable bead-spring models [2^ , where the migration velocity decreases with increasing 
rigidity of the tumbling object, corresponding also to increasing values of the spring constant. Indeed the vesicle 
deformability is besides the nonlinear shear gradient the main ingredient for the lateral migration in the Stokes limit. 
A vesicle in unbounded Poiseuille flow undergoes large deformations {Ca ^ 1, see Fig. [3]) caused only and mainly by 
the curvature of the velocity proflle. 



V. CONCLUSIONS 



The dynamical behavior of a single vesicle placed in an unbounded plane Poiseuille flow has been investigated 
numerically. We found that the vesicle migrates during its tank-treading motion towards the center of a parabolic 
flow proflle. The migration velocity is found to increase with the local capillary number (deflned by the time scale of 
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FIG. 4: (Color) A schematic showing the vector position r, the normal n and the tangent t unit vectors. 



the vesicle relaxation towards its equilibrium shape times the local shear rate) , but reaches a plateau above a certain 
value of the capillary number. This plateau value increases with the curvature of the parabolic flow profile c. When 
the vesicle reaches this final equilibrium position, its lateral migration velocity vanishes and it continues to move with 
a parachute shape parallel to the flow direction. We found that the migration velocity normalized to the curvature 
^^migration/c foUows essentially a universal law where the universal function depends on the local capillary number Ca, 
namely Wmigration/c ^ fiCa)- A droplet having no viscosity contrast seems to move away from the center which 
is in a marked contrast to vesicles, for which we found migration towards the center. This difference is not yet fully 
understood, and is currently under investigation. In forthcoming publications our calculations will be also extended 
to droplets in order to extract the main source of difference between the vesicle and drop migration. 
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In a two spatial dimension the vesicle membrane is represented by a one-dimensional closed contour. The corre- 
sponding membrane energy is an integral over this contour. 
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APPENDIX: DERIVATION OF THE MEMBRANE FORCE 




(A.l) 



where L is the vesicle perimeter (i.e. the length of the contour) and r the membrane vector position. Let 




(A.2) 



and, 




(A.3) 



The counterclockwise tangent unit vector (see Fig. [4|) is given by, 



t 

ds 



(A.4) 



and its derivative with respect to s defines the curvature. 




(A.5) 
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where n is the outward unit vector normal to the curve. The derivative of n with respect to s gives, 

'^=Ht. (A.6) 
OS 

Using Eq. (|A.4p and Eq. (|A.5|) we get the expression of the curvature, 



The membrane force is deduced from the functional derivative of the energy 5E/5v, where 5y is a local small 
displacement of the vesicle membrane. Due to the displacement of r by (5r, ds will undergo variations as well. It is 
convenient to introduce a fixed parametrization (instead of s) of the curve, which is denoted by a. a is a parameter 
that we can take to vary from to 1. The correspondence with s is such that s(q; = 0) = and s{a = 1) = L. We 
then introduce the metric g = \dr/da\ , so that ds ~ ^da. We convert the various terms in the energy by using 
now the variable a. The curvature assumes the following expression 



Writing d^v/da^ in terms of the tangent and the normal vectors, it is straightforward to show that, 

d^Y s 

This allows to eliminate s from the expression for H: 




1. The curvature force 

Replacing in Eq. (|A.2I by the expression given in Eq. (|A.11[) we obtain. 



Ec = ^l^' (^'~-^{^.if^g-'/'da. (A.12) 



The functional derivative of Ec reads (from classical variation results) 

SEc dec d dec dec 



(A.13) 



Sr dr da dr da^ dr ' 

whith ec = (k/2) ^r^ — ^ (fi")^^ 9~'^^^, i" — dv/da and r = d'^vjdc? . Since ec does not exphcitly depend on r, the 
first term on the right hand side of Eq. (|A.13p vanishes. The second term gives. 



d dec d ( 1 ( 1-:^^- , 3 ^..,2- 5 ,....2 



«7Jr + 7T(r)'r - TT-i^rYr , (A.14) 



da dr da \ g^/^ V 2g 

d f d^s H 3 
da \ do? g 2 



d ( d^sH 3. , 
^l^i-l^-^+7;Hh (A.15) 



while the third one becomes. 



d^ dec d^ f 1 .. 1 , 



da'^ dr da^ \g'^^^ g^l^ 

/ djHn) ^d^H\ 
da \ ds da^ g ' 
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Reporting the above results into (|A.13|) . we obtain the following expression for the functional derivative 

6Ec d f dH 1 



6r - ^a^^-a7"+2^ V' ^^-''^ 
Therefore, the membrane curvature force is given by, 

where the factor ^/g disappears from the physical force, since this one must be defined as fc = —{l/^)SEc/Sr, as 
explained at the end of the appendix. 



The tension force 



Finally Eq. (jA.3|) takes the following form 



whose functional derivative is, 



6Et 



d dei 



Sr da dr 

with er — C(i")V5- Note that er depends neither on r nor on r. We easily find 



SEj 



Sr 



da 



V9J 



-V^§; (C(r)t), 



-V9 



OS 



QHn 



The membrane force associated with the Lagrange multiplier is then. 



OS 



By adding Eqs. (|A.21|) and (|A.28p . we obtain the total membrane force 



f = 



9s2 



dC 



(A.22) 
(A.23) 

(A.24) 

(A.25) 
(A.26) 
(A.27) 

(A.28) 
(A.29) 



Let us briefly explain why the force is given by — f = —{1/^)5Et/5y (and not just —SEr/Sr). The reason is that 
what matters is a physical displacement of the curve element ds and not da (which is a mathematical arbitrary 
parametrization). If one performs directly the variation on the integral, one finds (according to the previous results) 



SE 



-HC 



n 



ds 
ds^ 



y/gdadr, 



dsSr = {dsSr 



(A.30) 
(A.31) 
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